Theory of atomistic simulation of spin-transfer torque in nanomagnets 



o 

(N 



i 

C 

o 
o 



> 
in 
t> 

cn 
o 



x 



Tiamhock Tay and L. J. Sham 
Department of Physics, University of California, San Diego, CA 92093 
(Dated: March 8, 2013) 

In spin-transfer torque (STT) for technological applications, the miniaturization of the magnet 
may reach the stage of requiring a fully quantum-mechanical treatment. We present an STT theory 
which uses the quantum macrospin ground and excited (magnon) states of the nanomagnet. This 
allows for energy and angular momentum exchanges between the current electron and the nano- 
magnet. We develop a method of magnetization dynamics simulation which captures the heating 
effect on the magnet by the spin-polarized current and the temperature-dependence in STT. We 
also discuss the magnetostatics effect on magnon scattering for ferromagnetic relaxation in a thin 
film. Our work demonstrates a realistic step towards simulation of quantum spin-transfer torque 
physics in nano-scale magnets. 

PACS numbers: 72.25.Mk, 75.75.Jn, 75.78.-n 



I. INTRODUCTION 

In an early study of the angular momentum transfer 
between a spin-polarized electron current and a ferro- 
magnetic thin film, 1 Berger modeled the fcrromagnet as 
a rigid classical spin which effectively acts as a spin- 
splitting field for the electrons. The torque acting on 
the magnet is then obtained by appealing to the classical 
angular momentum conservation law. 2,3 Berger acknowl- 
edged that this semi-classical treatment fails to capture 
the transverse quantum fluctuations. A recent study by 
Wang and Sham demonstrated how a fully quantum- 
mechanical treatment of spin-transfer torque is possible 
by considering an exchange interaction between each elec- 
tron and the magnetization, 4 the latter being modeled by 
a macrospin. Indeed, their study shows that restoring the 
transverse correlations between the magnet spins and the 
current spins gives rise to a noise in the magnetization 
dynamics of the magnet. This noise may become more 
important as the recording bits in magnetic storage me- 
dia approach the nano-scale regime. 

The goal of this work is to develop a theory of spin- 
transfer torque (STT) dynamics and a method of simula- 
tion suited to magnets down to nanometer size where the 
classical treatment of the magnet may fail. 1-3 We adopt 
the all-quantum approach of Wang and Sham 4 in treating 
both the current electrons and the magnet spins as quan- 
tum objects and, in addition, extend the treatment of the 
magnetization states beyond the rigid macro-spin states 
to include their excited states as magnons. The inter- 
action between the current electrons and the macrospin 
ground states involves only elastic scatterings 1-3 whereas 
the inclusion of magons accounts for energy transfer 
which is relevant to the thermal effects. 5 Our integra- 
tion of the magnon dynamics with the macrospin dy- 
namics includes the correlation between the transverse 
magnet motion and the current electrons, absent in the 
common magnon treatments in spin torque transfer 1 and 
the applications in magnonics. 6 We have also extended 
the primitive model of a delta function potential for the 
ferromagnet 4 to a slab to provide a simple model of dy- 



namics inside the ferromagnet for properties such as the 
thickness dependence. Thus, our theory provides a basis 
for an atomistic simulation of the quantum spin-transfer 
torque process for nano-magnets and the dissipation ef- 
fects of the magnons on the magnetization dynamics. 

An outline of the elements of our theory and key re- 
sults is as follows. In Sec. II, the Wang-Sham macrospin 
approach is applied to a model with a finite thickness for 
the ferromagnet to furnish the limiting case of our general 
theory where magnons are neglected. Sec. Ill contains 
an introduction of the magnon to the full Hamiltonian 
model for spin current driven magnetization dynamics 
and a description of the application of distorted-wave 
Born approximation to the current electron scattering 
with the localized spins of the magnet. The addition of 
the magnons in the scattering leads to energy as well as 
angular momentum exchanges between the current and 
the magnet. Sec. IV contains the formulation of a the- 
ory of wide-angle magnetization precession in an exter- 
nal field. Our theory of the magnetization-magnon in- 
teraction resolves the lack of relaxation mechanism in an 
insulating ferromagnetic film. 7-9 Using Monte Carlo sim- 
ulations to study the wide-angle dynamics, we find that 
magnetic impurity-induced damping results in deviation 
from the Landau-Lifshitz-Gilbert dynamics. Sec. V puts 
the nature of our theory and results in the perspective of 
the spin-torque field, especially in regard to future work 
along the concept of this paper. 



II. MACROSPIN MODEL WITHOUT 
MAGNONS 

We first establish the limiting cases of Refs. 1-4 by 
using Wang and Sham's approach in the coherent scat- 
tering of an itinerant electron by a thin ferromagnetic 
junction in one-dimension. The results of the scattering 
will be shown to reduce to the classical magnet case when 
the recoil terms are excluded. This work gives a clearer 
perspective of Wang and Sham's theory and serves as 
a limiting case test of our full treatment of spin torque 
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including spin wave excitations due to scattering by the 
itinerant electrons. Instead of a spin-spin interaction in a 
(5-plane, we use the following Heiscnbcrg exchange model: 



H s d — H e l + H int 
H el = I d 3 r ct 



ft 2 V 2 
2m 



-A (r) 



Hi nt = - J d 3 r A(r) j(r) • cr a pcl a c T p, 

~ -Jjf d3r A ( r ) J ' ^apC^Crp, 



(1) 

(2) 
(3) 
(4) 



where summation over each repeated index is implied, 
c ra are electron operators with spin component a, cr a p 
are Pauli matrices, and A(r) = A and Ao(r) = Ao in the 
ferromagnet but zero elsewhere. The magnet is described 
by a continuum background of physically reasonable spin 
density. Thus, j(r) in Eq. (3) represents a coarse-grained 
spin- 1/2 operator in the magnet. By replacing iVj(r) 
with the total spin operator J, where N = 2 J is the total 
number of localized moments in a fully-saturated magnet, 
Eq. (4) gives a macrospin model without spinwaves. In 
the limit of a macroscopic spin quantum number J, J 
may be replaced by a classical vector and Eq. (4) would 
then describe the interaction of an electron spin with an 
effective magnetic field A in the ferromagnetic film. 

The scattering can be block-diagonalized in the total 
angular momentum basis of the electron spin and the 
macrospin: 



S ■ J 



+J/2, F = J + 1/2 
-J/2, F = J-l/2, 



(5) 



where S is the electron spin operator and F is the quan- 
tum number of the total angular momentum operator 
F = J + S. Thus, it is more convenient to analyze the 
scattering in the respective total-spin sector 



H 2 V 2 
2m 



M --A(r)^(r), (6) 



where a = ± labels an electron wavefunction in the re- 
spective total angular momentum state with F = J+a/2. 
The interacting spin problem now reduces to a single- 
particle scattering problem. For convenience, we define 
the Berger frame as one in which the z-axis (Berger axis) 
lies along the initial magnetization direction of the mag- 
net prior to scattering by the electron. Consider a given 
input state in which an electron, with its spin pointing 
in the direction (6, 0), moves to the right toward the fer- 
romagnetic film. The scattering problem can be solved 
by first using Clebsch-Gordon coefficients to transform 
the spin basis to the total spin basis. After obtaining a 
scattering solution in the latter basis, we transform back 
to the former basis. 

For the following input state, 



a=±l/2 



| J, J) 



(7) 



we obtain an output state to the leading order in \j\fj 
expansion, 

|*out) ~ E xJ-(k',k)\k';l)®\J,J) 

k'=±k 



+ E l fc ';t)®|; 

fe'=±fc I 



xv 



Xt f+(k',k)\J,J) 

U(k',k)-f.(k',k) 



2 J 



\J,J-1) , (8) 



where f a (k',k) is the scattering matrix element (not 
scattering amplitude) for an electron in the total spin 
F = J + a/2 to scatter from momentum k to k' , cf. 
Eq. (15). This output state entangles the electron spa- 
tial and spin states, and the macrospin states of the mag- 
net. If two successive current electrons are uncorrelated 
prior to the encounter with the magnet, then the first 
scattered electron will decohere into a population dis- 
tribution of specific electron momentum and spin states 
which presents an associated mixed state of the magnet 
for the scattering of the second electron. This introduces 
noise in the magnetization dynamics of the magnet. 4 

In the infinite- J limit, the term drops off and 

Eq. (8) reduces to Berger's solution for an electron scat- 
tering off a rigid magnet. This is not unexpected since 
the Hamiltonian in Eq. (6) is formally identical to the 
scmiclassical counterpart if one replaces a = ± by fj,. In 
the classical approximation of the magnetization, 1,2 there 
is no recoil motion in the scattering and the recoil has to 
be inferred from the conservation of total angular mo- 
mentum of the current electron and the magnetization. 
The important consequence of the quantum treatment is 
that the recoil motion of the magnetization is present in 
the scattering terms of order 1/vJ. By comparing the 
second term with a spin coherent state [cf. Eq. (II)], 10,11 
the recoil in the direction of the magnet can be found for 
a given electron momentum and spin state. 

In Fig. 1, we use Monte Carlo simulations to study the 
effects of varying J and film thickness on the macrospin 
recoil. Using a model with a background spin density 
of 100 nm -3 , we compute for different film cross sec- 
tional area, the ^-component of the macrospin recoil as 
a function of the film thickness after a single electron 
(initially oriented in the — x direction) scatters from the 
magnet. For simplicity, we replace the discrete- valued J 
by a continuum J > 0.5 in this calculation. The recoil 
is averaged over all outcomes in Eq. (8). At w ~ 2 nm 
where /+ <C /- , we observe an appreciable difference be- 
tween Yong-Sham's macrospin recoil and that of Berger 
when J <C 10 3 . Thus, semiclassical treatments including 
micromagnetics remain valid for the mean magnetization 
dynamics down to the middle of the mesoscopic regime 12 
but are not able to account for the noise of quantum 
origin. 4 The thickness dependence shows that our quan- 
tum mechanical electron-magnet scattering is capable of 
capturing the band structure information in a multilayer 
device. 
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FIG. 1: [Color online] For a magnetic film with a given cross 
section, the rc-component of the macrospin recoil oscillates as 
a function of the film thickness when a single electron (ini- 
tially oriented in the —x direction) scatters from the mag- 
net. The recoil shown in the plot has been averaged over all 
the outcomes in Eq. (8). The result for J/w > 10 3 nm _1 is 
visually indistinguishable from that obtained using Berger's 
treatment. At w ~ 2 nm, quantum-classical differences are 
noticeable at the mesoscopic range J ~ 200 and prominent at 
the microsopic range of J ~ 20. 



III. SPIN-TRANSFER TORQUE WITH 
MAGNONS 



the spin and boson operators used to depict the local 
spins will be labeled with suffices for their positions or 
lattice momenta. Without loss of generality in represent- 
ing the dynamics of the magnetization, we may choose 
its instantaneous axis to make the HP expansion. Ref- 
erences 10,11 show that a macrospin state (Dicke state) 
may be expanded in terms of the coherent states and vice 
versa. Then, for J> 1, a correspondence between a bo- 
son coherent state and the macrospin state is established, 

|?;) =e-l"l 2 / 2 e" 6t - I '* b |0), (10) 

|J;M)=e- ,J * T \J,J), (11) 

(1 + ICI 2 ) 

i 1 = (V2J 1 C = e # tan(0/2). (12) 

The large J requirement is easily satisfied even for a small 
nanomagnet with 10 5 spin- 1/2 moments. For 2 J such 
spins, there are as many macrospin states, | J, J — l,r), 
r = 1, . . . , 2 J. For a lattice of local spins 1/2, Fourier 
transforms of these states give the spin waves or magnons 
&J|0), 14 with wave vectors q. The \J, J — l,q = 0) = 

6^ =0 |0) state is the state which provides the change in 
the rigid rotation of the macrospin state | J, J) through 
a small angle 9. This state will be used to determine 
the direction of recoil for a nanomagnet after an electron 
scattering. 



For a realistic treatment of the interaction between an 
itinerant electron and a nanomagnet, we add to the rigid 
macro-spin states of the magnet excited states approx- 
imated by magnons. In this section, we will first state 
the Holstein-Primakoff (HP) transformation and then re- 
late the correspondence between a rotated ferromagnetic 
ground state and a bosonic coherent state in the HP rep- 
resentation. We then introduce the Hamiltonian for the 
electron-ferromagnet system. By analyzing the interac- 
tion terms in the HP model, we identify the term corre- 
sponding to Berger's scattering potential and terms that 
provide quantum corrections to the electron-magnet scat- 
tering solution. In the distorted-wave Born's approxima- 
tion, we show that the solution of the preceding section 
is exactly reproduced by the HP model when spin wave 
excitations are neglected. The Born approximation is jus- 
tified by the small change in the excited magnetic state 
caused by one itinerant electron to the order of 1 / VJ. 

A. Spin coherent state in HP representation 

The Holstein-Primakoff representation for a spin-j op- 
erator is given by 13 

j-=bW2j-Vb, j+=jl, J, = J - 6+6, (9) 

where the spin state \j, j) serves as the vacuum state 
|0) for the boson operators b and . In the following, 



B. The full electron- magnet Hamiltonian 

We now restore magnons to the electron-magnet model 
by using the full interaction term Eq. (3) and also adding 
magnon energy to the Hamiltonian. Using the HP repre- 
sentation from Eq. (9), we make a series expansion up to 
quadratic order in the bosonic operators and obtain the 
following: 



H = H + V, (13) 




where w q = J'a 2 q 2 /2 is taken as a phenomenological 
magnon dispersion for the ferromagnet, with J a 2 be- 
ing the spin stiffness and a the spacing between the 
spin-1/2 moments which we have coarse-grained. For 
a given magnet with dimensions L Xl L y and L z , the 
magnon wave vector has components taking discrete val- 
ues qi = 2ixn.il Li where rij = 0, 1, . . . , Li/ a. We have 
neglected the magnon-magnon interaction since the elec- 
tron scattering time is much shorter than the relaxation 
time of a magnon mode. 
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(a) (b) (c) 

FIG. 2: Illustration of the interactions between an electron 
and a nanomagnet in the HP representation, (a) Classical 
scattering potential in Berger's treatment, (b) Electron spin 
flip causes orientation recoil of the magnet or spin wave excita- 
tions, (c) Electron-magnon interaction results in an effective 
renormalization of the classical scattering potential, leading 
to a temperature dependence in STT. 



q^O, a magnon is created. This leads to a heating of the 
ferromagnetic state by the spin-polarized current. The 
stochastic selection by a subsequent electron of the mag- 
netization states from the elastic and inelastic scattering 
by the preceding electron determines the new direction 
and the internal state of the magnetization. 

The second electron-magnon term in V further con- 
tributes a quantum correction of the following form 
|±k'; a) <£> 6q,6q \M). For q' = q, the magnon population 
has the effect of renormalizing the S z -j z spin coupling be- 
tween the electron and the magnet and thus contributes 
to temperature dependence in the orientation recoil. In 
the rest of the section, we develop a scattering theory to 
study the effects of the electron-magnet interaction. 



The spin-dependent term in Hq comes from the 
electron-magnet interaction; it gives the classical scatter- 
ing potential in Berger's theory and may be represented 
using the Fcynman diagram in Fig. 2(a). The effect of 
this term is discussed in Sec. II, where the recoiless scat- 
tered state is given by the leading-order term in Eq. (8). 

The two terms in V are depicted by Figs. 2(b) and 2(c). 
The cubic term contributes to Berger's scattering solu- 
tion a quantum correction of the form |±k'; f) <E> \M), 
where \M) denotes the initial state of the magnet prior 
to scattering. For q = 0, &t |0) = | J, J — 1) results in a 
recoil in the orientation of the magnet, and later we will 
show that this reproduces Wang and Sham's macrospin 
result exactly in the leading 1/yfJ approximation. For 



C. Classical scattering potential 

This subsection considers the one-dimensional problem 
of an electron scattering off a spin-dependent classical 
potential (i.e. without the interaction term V) which 
extends from x = —w/2 to x = w/2. Scattering solutions 
that will be used for calculations in later sections are 
those prepared with (i) an incident electron approaching 
the magnet from the left, (ii) a reflected electron moving 
left from the magnet, and (iii) a transmitted electron 
moving right from the magnet. It is a simple exercise to 
obtain the following solutions: 




k&x j ^ g ik(j 
ikx 



fa{+k, k)e 



+ f*(-k,k)e lkx 
g~e + h* a e lkaX 

-ikx 




g-^e""<*~ + h* a e lk " x 
„ikx + f*(-k,k)e~ ikx 



,U(+k, k) = -Akk a /c a , g a 
f a (—k, k) = 2i(k 2 - k 2 ) sm(k a w)/c a , h a 
h 2 kl/2m^h 2 k 2 /2m + S^ + (7X/2, c<7 



— oo < x < —w/2 
-w/2 < x < +w/2 
+w/2 < x < +oo 

— oo < x < —w/2 
-w/2 < x < +w/2 
+w/2 < x < +oo 

— oo < x < —w/2 
-w/2 < x < +w/2 
+w/2 < x < +oo 



(15) 



-2k{k + k a )e l(k - k " )w/2 /c a , 
+2k{k-k a )e l(k+K)w/2 / C(yi 
(k - k a ) 2 e l(k+k " )w - (k + k a ) 2 e l{k ~ k ° )w . 



(16) 



For the state $+k*( x ) in Ec h ( 15 )> .M±M) 

give the 

matrix element for a spin-er electron to be transmitted or 
reflected from the magnet. ^_k„(x) is the complex con- 
jugate of ^^(x), while \x) is the mirror reflection 

of ^_ fe (a;). In the following sections, these states which 
are commonly known as the "distorted waves" , will serve 



as zeroth-ordcr wave functions for the first Born approx- 
imation of the finite q magnon contributions. 
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D. Energy and angular momentum exchanges 

1. Distorted-wave Born approximation 

Although the exchange coupling constant A which ap- 
pears in Hq as well as V is not necessarily small, we 
explain below that it is still reasonable to treat V as 
a perturbation on Hq. Averaging over all N localized 
moments, (b\b r ) ~ n tot /N shows that c]. a c Ta b].b r indeed 
gives a small contribution when the temperature is much 
lower than the Curie temperature (Tc), since the total 
magnon population n tot <C N. For cL c r ^b r , only a small 
number of magnon modes with non-zero occupation con- 
tribute to the scattering when T <C Tc- The contribution 
from c\.^c r ]b\ is comparable to A, but for T <C Tp (the 
Fermi temperature), a large number of scattering out- 
comes with magnon creation are forbidden by Pauli ex- 
clusion principle since the outgoing electron has to scatter 
into a state deep below the Fermi level. 

Hence, we may treat V as a small perturbation. Ap- 
proximate scattering solutions can be obtained from the 
following T-matrix formula: 15 



T, 



an 



T 



(0) 



0a 



(4- } \vK +) ), 



(17) 



where a and (3 label the eigenstates of Hq, and and 
Tp a are T-matrix elements for an incoming state a to 
scatter into outgoing state /3 for respective Hamiltoni- 
ans Hq and Hq + V. Here, 4'q ± ' ) are exact scattering 
states of the Hamiltonian Hq, with +(— ) indicating in- 
coming (outgoing) states. Note that the total energy of 
the system has to be conserved by including energy ex- 
change between the electron and the magnon states. The 
first Born's approximation is used to derive the above 
formula. 15 



2. Recoil in orientation of the magnet 

We now determine the change in orientation of the 
nanomagnet, the magnon creation or annihilation ampli- 
tudes for inelastic events after electron dccoherencc oc- 
curs. For the electron-nanomagnet scattering model, we 
make use of exact expressions for TjfJ and \^^) from 
Eq. (15). Using the T-matrix formula, we derive explicit 
expressions as the starting point to investigate the effect 
of magnons on the rigid rotation of the magnetization. 
Here, we retain only the 6 q= o term in V while leaving 
the q ^ terms to the next section. The Hamiltonian 
relevant to the scattering between an itinerant electron 
and the macrospin is 



H = Hq 



A 



/ dx (6+ct Cri + h.c. , (18) 

N J-w/2 V ' 



where the number of spins N = 2 J. Taking overlap for 
the spatial part using the scattering solutions in Eq. (15), 



we obtain the following matrix elements: 

pw/2 



VJV J-w/2 

(19) 

_ n 2 fc/ T (±fc,fc)-/ i (±fc 1 fc) fet 



N 

For an incident electron with spin state (xt, Xl)> the scat- 
tered state of the electron-magnet system is 



E Xlh(k',k)\k';l)®\M) 



fc'=±fc 



+ E W;t)®\xtMk',k)\M) 

k>=±k { 



2 J 



-bl \M) 



(21) 



If the magnet is initially in the ferromagnetic ground 
state where \M) = \J,J) and b \M) = \J,J-1), the 
HP scattering solution reproduces Eq. (8) exactly. Thus, 
the distorted-wave Born approximation indeed allows an 
accurate treatment of the macrospin recoil in the leading 
1/VN order. 



3. Effect of magnons on orientation recoil 

Consider first the two-magnon interaction with the 
itinerant electron in V , 

^ E b lK+^ E ff / A A ( r ) e?q2 ' r c -c rCT . (22) 
qiq2 <y 

For q2 = 0, the contribution rcnormalizcs the S z -j z cou- 
pling between an itinerant electron and the magnet. Ne- 
glecting terms with q2 ^ 0, we move the diagonal terms 
to the non-interacting part of the Hamiltonian Hq and 
replace the magnon occupations by their mean-field ex- 
pectation values. This approximation leads to a renor- 
malized, anisotropic spin coupling between the electron 
and the macrospin: 

ci a {cr a p ■ j r )c rf ) ~ 4 t c ri ,j~ + 4 4 c rt j+ 



a -cl a c ra |l-^W/jJ . (23) 



The scattered state in Eq. (21) is accordingly extended 
to 

|*o«t) « J2xJl(±k,k)\±k;l)®\M) 
± 

+ El± fc ;t)® jxt/t(±M)|M) 



Xi 



M±k,k)-h(±k )k ) blm \ m 
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time (arbitrary unit) 

FIG. 3: At high temperatures, the reduced total spin of the 
magnet should lead to a larger orientation recoil as electrons 
scatter off the magnet. But a smaller difference in the scat- 
tering amplitudes for spin up and spin down electrons at high 
magnon population over-compensates and results in a reduced 
STT. The inset shows a vertical cut at t = 0.32, plotted 
against the magnet temperature. The scatter results from 
the stochastic nature of individual Monte Carlo runs. 



where f a (±k,k) are the matrix elements of an electron 
scattering from a magnet, computed using the reduced 
spin splitting A[l — X) q ( n q)/^1- Two factors with com- 
peting effects are present in the spin flip term. The ex- 
tra factor [1 — ^2q( n q) / J] hi the denominator acts to 
increase the recoil magnitude at large magnon popula- 
tions, while the reduced S z -j z coupling leads to a smaller 
difference between the spin up and spin down scatter- 
ing amplitudes. To determine the resulting temperature 
dependence in magnetization switching, we simulate the 
trajectories of the magnetization at various temperatures 
using Monte Carlo. Figure 3 shows that the reversal 
rate decreases with increasing temperature, thus allow- 
ing us to conclude that STT is weaker in the presence 
of magnons. Note that this does not conflict with heat- 
assisted magnetic recording (HAMR) since the latter oc- 
curs via a temporary loss of magnetism by heating a mag- 
net beyond its Curie point, followed by remagnetization 
in a magnetic field as it cools. Note also that we have ex- 
cluded magnetic and anisotropy fields in this simulation 
so that the activation effect of an energy barrier on the 
switching is removed and hence, allowing us to isolate 
the effect of magnon population on the STT. Therefore, 
Fig. 3 should not be confused with "thermally-activated 
switching" ; the latter pertains to the Arrhenius factor 
in magnetization reversals. 16-18 Full treatment of STT- 
driven switching in a field will be deferred to Sec. IV. 

We now discuss the neglect of spinwavc transforma- 
tion under rotation in our formulation. In Eq. (24), one 
observes that the operator 6q rotates the ground state of 
the magnet but not the spinwaves. For small-angle recoils 
typical of nanomagncts with at least 10 5 localized spin- 
1/2 moments, an 0(N~ 1 ^ 2 ) correction to the spinwavc 



operators under rotation is indeed negligible. Therefore, 
we may regard b Q \ M) in Eq. (24) as a rigid rotation of 
all the spins in the magnet. Alternatively, one may ap- 
proximate the state of a magnet at finite temperature by 
| J', J') with a reduced total spin J' = J[l — X) q ( n q)/^] 
and analyze the orientation recoil by appealing to the 
original interaction Hamiltonian: 




(26) 

The first term in Eq. (26) causes a rigid rotation of the 
spinwaves along with the uniform magnetization as an 
electron scatters off elastically from the magnet. The 
remaining terms involve inelastic processes and will be 
discussed next. In Born's approximation, the orientation 
recoil is analytically identical to the preceding formula- 
tion at the leading 1/N order. Thus, the HP represen- 
tation allows one to study magnetization recoils at zero 
temperature as well as finite temperatures. 



4- Magnon creation or annihilation 

The remaining single-magnon and electron interaction 
terms for the derivation of the coefficients of amplitudes 
of states with a magnon created or annihilated, 19,20 are 




As before, we use the scattering solutions in Eq. (15) to 
obtain the following matrix elements 

(^imff) = -^E' w n(k',k;q)6 q , (28) 

= --L5>it(k',k;q)& q , (29) 
* q 

<w(k',k;q) = / dx (x)e*«"*,(30) 

J-w/2 

where the primed sum indicates restriction of the scat- 
terings to an equal total energy shell. Then the magnon 
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500 




time (arbitrary unit) 

FIG. 4: Monte Carlo simulation of the spin lattice tempera- 
ture during magnetization switching. Due to a large number 
of magnon modes with zero occupation number, the proba- 
bility for magnon creation is greater than that for absorption, 
thus leading to an overall monotonic temperature rise during 
the switching process. 



generally higher for magnon creation. For relative orien- 
tation of the magnetization and electron spin close to the 
anti-parallel case, \xi\ 3> |xtl results in more magnon 
creation than magnon destruction. When the electron 
spin and the magnetization are nearly parallel, the oppo- 
site occurs. However, the destruction of many low-energy 
magnons together with the creation of fewer high-energy 
magnons do not necessarily result in a net heat extraction 
from the magnon bath. This explains the monotonous 
temperature rise of the magnon bath. On average, less 
than 5% of the scatterings in our simulation are inelas- 
tic; this justifies the distorted-wave Born approximation 
used. For this simulation, we exclude couplings between 
the magnon, electronic and phonon degrees of freedom. 21 
Inclusion of these effects phenomenologically should re- 
sult in a temperature peak followed by a cooling of the 
spin lattice during the magnetization reversal. 



IV. FIELD EFFECTS ON MAGNETIZATION 
DYNAMICS 



correction to the scattering solution is 

q,k 

+ F ^y^t(k',k;q)|k / ;4)®6 q |M)|, (31) 

k'(E,q) = ^^T^mujjh 2 . (32) 

Equation (31) shows that a spin-down electron flips up- 
ward by emitting a magnon while a spin-up electron flips 
downward by absorbing a magnon, with the energy bal- 
ance provided by the kinetic energy of the electron in 
Eq. (32). 

Figure 4 is a simulation of the time evolution of the 
magnet temperature as a nanomagnet undergoes STT- 
driven switching. After each inelastic scattering event 
that occurred during our Monte Carlo simulation, we fit 
the non-equilibrium distribution of the magnet states by 
a single effective temperature before the next electron is 
injected into the magnet. More rigorous treatment of 
the magnetization dynamics will be left for future stud- 
ies. The result shows that the spin temperature of the 
nanomagnet rises monotonically on average, but with a 
decreasing rate of rise. This is consistent with the scat- 
tering result in Eq. (31). For magnon absorption, the 
contribution comes only from low energy states with non- 
zero occupation. On the other hand, all magnon states 
contribute to magnon creation regardless of their occu- 
pation numbers. Thus, the total probability for an elec- 
tron to emit a magnon is evidently much greater than 
that for absorbing a magnon. Inclusion of the effect of 
Pauli exclusion on an electron scattering into states below 
the Fermi level of the normal metal leads to comparable 
probabilities for magnon creation and destruction, but 



In the preceding sections, we developed a consistent 
theory of energy and angular momentum exchanges be- 
tween an itinerant electron and a magnet by treating the 
magnet states quantum mechanically. We now turn to 
the problem of field-driven precession in the intervening 
time interval between two successive electron scatterings. 
Consider a fixed reference frame that coincides with the 
Berger frame at time t = 0. At any time, the ground 
state of the magnet is always given by the vacuum state 
in the instantaneous Berger frame, but relative to the 
fixed reference frame, |0) transforms into a coherent state 
|jjo(i)> defined in Eq. (10). By Vo (t) = M*)l&o|»to(t)> 
and the Heisenberg equation for &o, we obtain fjo{t). Al- 
though the HP approach is restricted to small-angle dy- 
namics about the quantization axis, the full trajectory 
of the magnetization can be built up using infinitesimal 
time evolution in a series of Berger frames. For the pur- 
pose of modeling, one might approximate the quantum 
state of the magnet using a direct product of magnon 
coherent states. This avoids the need for a prohibitively 
large amount of resources in simulating the distribution 
of magnet states. We leave the simulation of the q ^ 
magnon evolution for future studies. 



A. HP approach for wide-angle precession 

In addition to the Zeeman energy of a magnet in an 
external magnetic field, we include the effect of dipolar 
interaction, crystalline and shape anisotropies by defin- 
ing an effective field h that exerts the same net torque 
on the magnet. 22 To treat relaxation in the magnetiza- 
tion dynamics, we include a b^bq scattering term in our 
Hamiltonian and obtain the spin relaxation in terms of 
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rates. The starting magnon Hamiltonian is 
W^b + - k 

q#0 



V „ - n V 



+ w q )& q fe q - ^— (h+bl + ft_&o), (33) 



where h± = h x ± w q contains contributions from 
Heisenberg exchange and demagnetization field, 8,13 and 
Wq gives the interaction between the uniform mode and 
the magnons. In Ref. 23, Sparks et al showed that the 
Wq term arises from impurity fields and is responsible for 
damping. The terms containing h± or h z are obtained 
by expanding the HP representation of the effective Zee- 
man energy —h - J to cubic order and then transforming 
to momentum space. We exclude the current electron- 
magnet interaction from Hm and apply this Hamiltonian 
in between two successive electron scatterings. In the 
following, rjq represents the magnon expectation in the 
coherent state |?7 q (f)). Taking expectation of the Heisen- 
berg equation of motion for bq with respect to |f/ q (£)), we 
obtain 



(34) 



V N 

{h z + w q )?7 q - 1 — h + 8qp 



2VN 



5>* 



?7q+q' 



q' 



q' 



Note that the constant 0{ v N) term results in a strongly 
driven harmonic oscillator coupled weakly to a reservoir 
of interacting oscillators. Ignoring all terms except this, 
we obtain dtrjo(t) = i\f~Nh + /2. Equation (10) gives 
Qe 1 ^ = ih + which results in a trajectory pointing along 
z x h in the Berger frame. In a fixed reference frame, this 
reproduces the classical result d^M = M X h. 

To obtain the damping effect due to the magnons, we 
use ??o = f lQ e-' l{ - h - +Q ^ t + Vn^q to obtain 



9e^ = 

r = 



i(2h z + 2Cj q - iT)£ , 



i 



W = Wo 



(1-ieoi 2 ) 3 / 2 

with ^o = \£o/h+\h-\. satisfying 



\*ioh + \, 



(36) 
(37) 

(38) 



2(1- |€oP)3/ 2 + \^\ - 1 - ^ U « /N - m 



The derivation of Eqs. (36)-(39) are given in Appendix A. 
Since the precession frequency shift due to anisotropics 



and demagnetization are already captured by the effec- 
tive field h z , 22 we set ujq to zero and obtain 



ih + (l-n to t/N) +£oI\ (40) 

0.5/h z , \h+\^+h z . 
0.6/|/i + |, \h+\*\h,\. (41) 
i.o/\h+\, \h+\<:-h,. 



Equation (40) shows that the frequency of precession 
decreases with the total magnon population n to t, and 
is therefore strongly dependent on temperature. 24 We 
have neglected the frequency detuning due to the Wq 
scattering term [see text after Eq. (A4)] and instead fo- 
cus on the damping effect that causes the magnetization 
to spiral toward the field. In a fixed reference frame, 
Eq. (40) reproduces the Landau-Lifshitz (LL) equation 

d t M = Mx(h-aMx h) for n tot = 0,witha = \£ T/h+\. 
Thus, the phenomenological damping parameter a in the 
LL equation is in fact not constant, but dependent on 
the angle 9^ between the magnetization and the field h. 
The LL equation grossly underestimates the true damp- 
ing when 9h is close to tt. Since a also depends on Wq, 
a microscopic derivation of Wq would be necessary for 
studying the relaxation physics beyond the LL theory. 
We will not go into such details here, 25 but instead ex- 
amine the magnetization dynamics resulting from specific 
field-dependence of a in the next section. 

The delta function in Eq. (37) indicates energy con- 
servation when exciting a magnon mode. The nonlinear 
frequency renormalization in Eq. (38) does not appear 
in the semiclassical treatment of the ferromagnetic reso- 
nances. For a Heisenberg ferromagnet in a magnetic field, 
the uniform mode lies at the minimum of the magnon 
dispersion curve and therefore, has no relaxation. An- 
derson and Suhl showed that magnetic dipole interaction 
causes the magnon dispersion curve to split into a band 
which depends on the magnitude and direction of the 
magnon wave vector, 7 thus creating some phase space 
for satisfying the constraint ujq = cj q . 8 This degenerate 
magnon mechanism, as it was known, does not work for a 
magnetic thin film when the magnetization has a signif- 
icant out-of-plane component. 9 Equations (38) and (39) 
show that an additional shift of 0(\h\ sin 2 9 h) or 0(\h\) is 
present respectively for 9h less than or greater than tt/2, 
where 9t is the angle between h and the Berger axis. This 
was not found earlier since the field direction is chosen 
as the z-axis in conventional HP expansion (9^ — 0) to 
eliminate the linear term in 6q. 



B. STT-driven dynamics in a field 

With a quantum approach for modeling STT-driven 
magnetization dynamics in an effective magnetic field, we 
describe the implementation details for a Monte Carlo 
simulation for the study of STT-driven precession and 
magnetization switching in a magnetic film. 
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1. Description of Monte Carlo simulation 

Consider the current as a stream of electrons incident 
on the magnetic film from one side at a uniform rate 
and constant momentum, which is a reasonable simpli- 
fication due to insensitivity of the noise to the electron 
momentum distribution found in Ref. 4. For each inci- 
dent electron, we compute the scattering amplitudes for 
each of four possible elastic outcomes as well as the nu- 
merous inelastic outcomes. Occupation numbers for the 
electron levels in the normal metals are sampled stochas- 
tically according to their Fermi distributions and forbid- 
den outcomes are then projected out as required by the 
Pauli exclusion principle. Next, we sample a scattering 
outcome and determine the new orientation of the mag- 
netization in the case of an elastic scattering, or update 
the internal energy of the magnet for an inelastic out- 
come. The magnetization orientation is then allowed to 
precess under the effective field until the next electron 
impinges on the film. We assume that the magnet is 
in quasi-equilibrium in between successive electron scat- 
terings, and the temperature of the magnet is updated 
by requiring that the new total energy of the magnet is 
equal to that given by the Bosc-Einstcin distribution of 
magnons in the magnet. 

In the following, we illustrate our approach by apply- 
ing it to a magnet in various fields. We consider a mag- 
net comprising a 100 x 100 x 10 spin- 1/2 moments with 
J = 80 meV and A = —0.7 eV. All incident electrons 
have the same momentum k x = 9 nm _1 and an average 
spin polarization of 25% with respect to a certain direc- 
tion. The coupling strength between an itinerant elec- 
tron and the spin-1/2 moments in the magnet is taken to 
be A = 1.2 eV. In principle, we should implement Lan- 
dauer's principle numerically by sampling electrons from 
both sides of the film according to their respective Fermi 
distributions. But this is computationally very costly 
and we leave that for future studies where the noise in a 
much smaller magnet is instead the main focus. 



1.0 



0.5 



0.0 



-0.5 - 



-1.0 



a LL = 0, r = 0.02h E 



0.00 



0.05 



0.10 



0.15 



FIG. 5: Jz/J versus current for an STT-driven precession in 
a vertical magnetic field with damping a — olll + \£oT/h+\. 
The LL damping (qll = 0.01, T = 0) results in an abrupt 
transition between no-switching and complete-switching at a 
critical current density 0.017 A/)j>m 2 , while the damping with 
ocll = and V = 0.02fts shows a gradual switching for cur- 
rent density above the critical value. 



Let us define the Berger frame such that its z-axis is 
specified by (Qm,4>m) with respect to a fixed reference 
frame, and its y-axis lies along y cos (f>M —i sin 0m in the 
fixed frame. Using Eq. (40) to obtain the new magne- 
tization direction (6, </>) in the initial Berger frame, the 
new orientation {9 / m -,4''m) m the fi xe( i frame is given by 

cos6' M = cos 8m cos 9 — cos sin 6* m sin 6*, (43) 



t> M + arg 



cos t) sin om 



+ sin 9 (cos 0m cos § + i sin (/>) . (44) 

Note that the Berry phase associated with the rotating 
Berger frame is implicitly captured by these equations. 



2. Effective field 

In the Berger frame, the effective field experienced by a 
single-domain magnet, due to the Zeeman energy — hp • J, 
an easy-axis anisotropy energy hx sin 2 9k, and an easy- 
plane anisotropy energy hp cos 2 8p, is 

h = Kb + h-K cos 9k — hp cos 9p, (42) 

where hK and hp respectively lie along the easy- axis and 
the normal to the easy-plane, 9k is the angle between hK 
and J, and 9p is the angle between hp and J. The de- 
magnetization field of a thin magnetic film can be lumped 
into hp. For dimensionlcss spin operator J, each field is 
given in units of energy. Thus, 1 kOe would correspond 
to 11.6 /xeV for the electron gyromagnetic ratio 5 = 2. 



3. Magnetic field only 

Consider an external field with hp = 1 kOe along the 
negative z-axis and a current with average spin polariza- 
tion along the z-direction. With this neglect of the de- 
magnetization field and all anisotropies, Fig. 5 compares 
the J z component of the processing magnetization for dif- 
ferent damping a, with a = aLL + \£oT/h+\. For Landau- 
Lifshitz damping with cull = 0.01 and Y = 0, com- 
plete switching occurs for any current above 0.017 A//im 2 
(Landau-Lifshitz-Gilbert equation would give similar re- 
sult since a is small). For cull = and T = 0.02hp (the 
magnitude of T is chosen to give the same critical cur- 
rent), J z switches progressively after the current crosses 
the critical value. 
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FIG. 6: Precession frequency against current through a mag- 
netic film (with electron spin polarization normal to the film), 
for damping a = oll + \£oY/h+\. The film contains easy- 
plane anisotropy but no other anisotropy or field is present. 
The precession frequency for «i damping (gill = 0.01, T = 0) 
is proportional to the current, while the ai damping (ull = 
0, r = 0.01/ip) imposes a threshold current for exciting a 
precession, and the 03 damping (qm = 0.008, Y = 0.002/ip) 
lies between the two. The precession frequency saturates at 
13.1 GHz. 



4- Easy-plane anisotropy only 

Here, we consider a magnetic film with an easy-plane 
anisotropy hp = 5 kOe (inclusive of a contribution from 
the demagnetization field), but no other anisotropy or 
field is present. When a current with spin polarization 
along the hard-axis passes through the film, an STT- 
driven precession occurs with a frequency that increases 
with the out-of-plane magnetization. Figure 6 compares 
the precession frequency for different damping a, with 
ot = cull + |£c>r//i+|. For the LL damping (ull = 0.01, 
r = 0) the precession frequency is proportional to the 
current, indicating that the out-of-plane magnetization 
is proportional to the current. This would be observed 
if the Gilbert damping from adiabatic spin pumping ef- 
fect of Ref. 26 is much greater than impurity-originated 
scattering between the uniform mode and the magnons. 
The damping with c*ll = and T = 0.01/ip imposes a 
threshold current for exciting a uniform precession and a 
large current is required for complete switching. In exper- 
iments, one might perhaps observe the third case where 
olll = 0.008 and Y = 0.002/ip are chosen for illustration. 



5. Thin film with in-plane uniaxial and magnetic fields 

Lastly, we consider a magnetic film with hp = 7.0 kOe 
along the z-axis, hn = 1.2 kOe along the x-axis, and 
an in-plane magnetic field hp = 2.0 kOe applied at 5° 
to the easy-axis. Using Y = 0.8 fieV, Fig. 7(a) shows 
a Monte Carlo simulation of the magnetization trajec- 
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FIG. 7: (a) Magnetization trajectory driven by a current den- 
sity jQ = 0.1 A/fim 2 . The orientation spirals away from its 
initial direction before reaching a steady orbit, (b) For j'q < 
0.078 A/^m 2 , the magnetization undergoes a small-amplitude 
precession at a constant frequency 15.2 GHz. Above this crit- 
ical current, the trajectory jumps to a large, distorted orbit 
with a precession frequency that decreases with the current. 

tory driven by a current density j'q = 0.1 A//im 2 . The 
strong easy-plane anisotropy results in a highly distorted 
trajectory. Figure 7(b) shows the dependence of the 
precession frequency on the current. Below the critical 
current density 0.078 A//mi 2 , the magnetization under- 
goes a small-amplitude precession with a constant fre- 
quency 15.2 GHz. Above this critical current, the tra- 
jectory jumps to a large, distorted orbit with a preces- 
sion frequency that decreases rapidly as the current in- 
creases further. This qualitatively reproduces the classi- 
cal Landau-Lifshitz-Gilbert macrospin simulation in the 
experimental study of Ref. 27. 



V. DISCUSSIONS AND CONCLUSION 

We have presented a simple microscopic model in a 
given field, together with contact exchange interaction 
between an itinerant electron and the spins in the mag- 
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net and detailed the quantum scattering between the 
current electron and the rigid macrospin states and the 
excited states (magnons) of the nanomagnet which pro- 
duces SST-driven precession and switching. The quan- 
tum approach in the large J limit produces the current 
electron spin rotation which is used with angular mo- 
mentum conservation to infer the macrospin recoil of 
Slonczewski 2 and reproduces the zero-wave vector spin 
wave emission of Bcrger. 1 By contrast, the quantum 
method deduces the recoil dynamics directly without the 
explicit use of angular momentum conservation. The use 
of angular momentum conservation for the spin torque 
transfer is valid where there is a large number of incoher- 
ent degrees of freedom involved but may be broken in ap- 
proaching the quantum limit through uncertainty created 
by either ultrafast control or atomic scale spin localiza- 
tion. Our simple simulations demonstrate the capability 
of the quantum approach to study damped magnetization 
dynamics starting from a microscopic treatment of ferro- 
magnetic relaxation, the noise sources from the trans- 
verse spin fluctuations in unitary scattering and energy 
transfer from the inelastic scattering involving magnons. 

We note that an important difference exists between 
our Holstein-Primakoff boson expansion and existing 
studies in the literature. The latter often expand the spin 
operators about a fixed magnetic field, thus restricting 
studies of magnon relaxation to small angles between the 
magnetization direction and the applied magnetic field. 
In constrast, our expansion in Eqs. (14) and (33) does 
not require the magnetization direction to be close to a 
fixed axis. We showed that the electron-macrospin inter- 
action in our HP expansion exactly reproduces Wang and 
Sham's quantum correction at the leading 1/y/N order. 
Furthermore, the HP model is able to treat precession 
damping even at large angles from the field. 

We have treated the current electrons as uncorre- 
cted with one another in the scattering events with the 
magnet which provide a stochastic source for the STT- 
driven magnetization dynamics. In our simulations, we 
left out a few aspects that are to be treated quantum- 
mechanically. First, fluctuations in the macrospin di- 
rection in connection to the fluctuation-dissipation the- 
orem for macrospin-magnon scatterings in a field had 
not been included. Second, the macrospin recoil due 
to adiabatic spin pumping effect of Ref. 28 was also 
excluded. Third, contributions from electron-magnon 
scatterings to relaxation in ferromagnetic metals had 
not been considered. Together, these can be taken 
into consideration by replacing the damping term with 
a phcnomenological Wiener process with drift, or one 
might try to extend our quantum approach to include 
them. All these would be necessary for simulating the 
"thermally-activated switching". 17 Our solution of the 
magnetization dynamics also has to be extended to ap- 
ply to more coherent sources of the driving mechanism, 
such as microwave, 27 laser, 21 and external spin waves. 29 
The quality of the quantum approach also has to be 
tested for nano-magncts or magnetic dots in the areas of 



spin pumping in metallic 27 or insulating fcrromagnets, 28 
magnonics 30 , and domain wall dynamics, 31 in which the 
magnetization had hithertofore be treated as classical us- 
ing micromagnetics. 
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Appendix A: Derivation of Eqs. (36)-(39) 

First, we substitute 770 = r/ + 770 into Eq. (35), with 
constant 770, and obtain 

q q 

+ (h z + u )fjo + -^j=(rfov'o + V0V0 + Vov'o* ) 
2v Jy 

+ ^(^^)^E'^, (AD 

v v q 

Note the grouping of terms linear in r)' that contribute 
to a renormalizcd energy for the uniform mode. The re- 
maining terms are either approximately constant or have 
phases that rotate with frequencies far from that of r]' . 
Since 770 = in the Berger frame at t = 0, we further 
substitute r] = —770 = —yNh+To (with real To) into the 
last two lines of Eq. (Al) and obtain 

idtVo = ^Nh + {\h + \ 2 rl + (h z +u )r - ± + 

q 

+ (^+ Wo + |/V| 2 T ) 77 + ^EV q *77 q . (A2) 

* q 

We have neglected 77 q 77_ q since this gives a higher order 
correction to damping compared to the Wq term. Wc 
now eliminate the strong driving term on the first line 
of Eq. (A2) by an appropriate choice of To. Defining 
r fo = 77oe~ i (' l2+ "°' t and 7y q = 77 q e _ "" qt , where Q = oj + 
\h + \ 2 r , we use the Weisskopf-Wigner approximation to 
obtain 32 

= -TjEVql 2 I dt'e^-^', (A3) 

iV q J 

« -VoT/2, (A4) 

where T is defined in Eq. (37). The approximation made 
in Eq. (A4) consists of extending the upper limit of the 
integral to infinity and neglecting the principal part of 



12 



/ °° dt e lLUt = iP(l/ui) + it5(uj). Transforming back to 770 
then gives Eq. (36), where £o in the main text is related 
to To via £o = h + TQ. 

When the angle Oh between the magnetization and the 
field is close to n, Eq. (A2) gives a divergence in £ , 
thus indicating that the truncated HP expansion at cubic 
order is no longer valid. To remedy this, we include every 
term in the infinite HP series that contains only 770 or t]q. 
This gives 



ldtV0 = ~m^ w ^ 



Nh A 



7t+l (2n-3)!!?7g"??g 
2™ n! N n 



iV7j_ g n (2n - 3)!! v* n ^Vo +1 



2 ^ 2 n n\ 

1 1— 1 



VqVq + (h z 



2VN 



N' 



(A5) 



' + 



2VN 



q 

^E V q>q> (A6) 



U) 



q 

00 

W „ + E 

n=l 



n(n+l) (2n-3)!! |£ | 2 ™ 



2" 



7" 



(A7) 



By differentiating \J\ — x and its series expansion, we 
obtain the following identity, 



^n(n+l)(2n-3)» xn _ x(4-3x) ^ 



^ 2 

n=l 



which then gives 



(l-x) 3 / 2 ' 



where n!! = n(n — 2) ••■I, I = 1,2 for odd and even 
n respectively, and nil = 1 for n < 0. As before, we 
substitute 770 = Vo + Vo arL d Vo = V~N£,o = VNh + To, and 
obtain 



Nh A 



(n + 1) (2n-3)!! 



I^ol 2 "} 



with a) and ^ satisfying Eqs. (38) and (39) respectively. 
One can show that Cjq > ojq holds strictly. 
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